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The vertical structure of freely evolving, continuously stratified, quasi-geostrophic flow 
is investigated. We predict the final state organization, and in particular its vertical 
structure, using statistical mechanics and these predictions are tested against numerical 
simulations. The key role played by conservation laws in each layer, including the fine- 
grained enstrophy, is discussed. In general, the conservation laws, and in particular that 
enstrophy is conserved layer-wise, prevent complete barotropization, i.e., the tendency 
to reach the gravest vertical mode. The peculiar role of the /3-effect, i.e. of the existence 
of planetary vorticity gradients, is discussed. In particular, it is shown that increasing /3 
increases the tendency toward barotropization through turbulent stirring. The effective- 
£>~>j ness of barotropization may be partly parameterized using the Rhines scale 2irE^ 4 / /3 1 / 2 . 

As this parameter decreases (beta increases) then barotropization can progress further, 
because the beta term provides enstrophy to each layer. However, if the beta effect is too 
large then the statistical mechanical predictions fail and wave dynamics prevent complete 
barotropization. 



1. Introduction 

The vertical structure of eddying flow in the oceanic mesoscale is a fundamental prob- 
lem in geophysical fluid dynamics, one that has has been reinvigorated by th e need to 
interp ret altimetric observation of surface velocity fields ( Scott fc Wane! 2005 : Lapevre 
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2009). Energy cascades and self-organization processes are fairly well understood and 
characterized in two-dimensional turbulence, but how these results generalize to rotat- 
ing, continuously stratified flows is still an open question. We ask in this paper what 
is the final organization of a freely evolving inviscid quasi-geostrophic flow with con- 
tinuous stratification. In particular, we examine what are the precise conditions for the 
oft-cited barotropization process to occur. Barotropizatio n refers to the tendency of a 
quasi-geostrophic flow to reach the gravest vertical mode (|CharnevHl97l lRhinesHl977l) . 



Because the gravest vertical mode is the barotropic one, i.e. a depth independent flow, 
"barotropization" means a tendency toward the formation of depth independent flows. We 
study in particular the key role played by the beta effect (the existence of planetary vor- 
ticity gradients) in such barotropization processes. It has been previously noticed that the 
pre sence of these la r ge sca le planetary vorticity gradients favors barotropic process, see 
e.g. ISmith fc Vallisl ( 2001 ). but this phenomenon remains to our knowledge unexplored 



and unexplained. 

Continuously stratified quasi-geostrophic flows take place in three dimensions, but 
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their dynamics are quasi-two-dimensional because the non-divergent advecting veloc- 
ity field has only horizontal components. The dynamics admits similar conservation 
laws as the two-dimensional Euler equations, including among others the total energy 
and the enstrophy of each layer. These conservation laws imply an inverse energy cas- 
cade toward large scales, and w eak energy dissipation in presence of weak viscosity 
( Kraichnan fc Montgomervfl 980). This situation contrasts with three dimensional turbu- 
lence where the energy cascades forward toward small scales, giving rise to finite energy 
dissipation no matter how weak the viscosity. 

Just as for two-dimensional Euler flows, inviscid, freely evolving continuously stratified 
quasi-geostrophic flows are disordered and an unstable initial condition rapidly evolves 
into a strongly disordered fine-grained flow with filaments being stretched, folded, and 
consequently cascading toward smaller and smaller scales. However, the inverse energy 
cascade leads to the formation of robust, large scale coherent structures filling the domain 
in which the flow takes place, contrasting strongly with the analogous three dimensional 
flow with its energy transfer to small scales. 

The aim of this paper is to provide a prediction for the vertical structure of the large 
scale flow organization resulting from inviscid, freely evolving continuously stratified 
quasi-geostrophic dynamics using a statistical mechanics approach. The goal of the sta- 
tistical mechanics approach is to drastically reduce the complexity of the problem of 
determining the large scale flow organization to the study of a few key parameters given 
by the dynamical invariants, as for instance the total energy and the fine-grained en - 
strophy in each layer . Such a th eory has been independently proposed bv lRoberd (|l990h : 
Robert fc Sommerial ( 199lh and lMillerl (|l990f h and henceforth will be referred to as MRS 



theory, ft is an equilibrium theory, developed in the case without forcing and dissipation. 
Our contribution will be to apply this theory to the explicit computation of a class of 
equilibrium states of the stratified quasi-geostrophic equations. 

A crucial ingredient of the theory comes from the constraints given by dynamical 
invariants. In the present paper, we will focus on the role played by the energy and the 
enstrophy constraints. In the framework of the MRS theory, it is justified in a low energy 
limit, and it is only a first step before more comprehensive studies including the ef fects 
of additional invariants, see e.g iMaida fc Wane! (j2006l) ; ICorvellec fc Bouchetl (|201lh . 

The phenomenology of two-dimensional and geostrophic tur bulence is often explained 
by con s idering the ene rgy and enstrophy constraints only, see e.g lKraichnan fc Montgomery 
( 198Cll ): Salmon ( 19981 ). Actually, earlier statistical mechanics approaches were developed 
for Galerkin-truncated Euler or quasi-geostrophic dynamics, which al lowed for a s imple r 
theoretical treatmen t than in the continuous dynamics case, see e.g. iKraichnan (1967); 
Salmon et al. ( 19761 ). In the case of Galerkin-truncated models, only the energy and the 



enstrophy are conserved quantities, so the statistical mechanics of these systems were 
called energy-enstrophy theory. 

A first reason to consider MRS theory rather than earlier energy-enstrophy statistical 
mechanics approach is that real flows are continuous, not Galerkin-truncated. Because 
real flows are also forced-dissipated, our first working hypothesis is that the computation 
of the MRS equilibrium states of the continuous dynamics is a necessary step before 
taking into account t he effect of weak forcing and dissipation. This is for instance the 
approach followed by iBouchet fc Simonned ((2009) , who found that in the presence of 
weak forcing and dissipation, the dynamics were close to MRS equilibrium states which 
were different from energy-enstrophy equilibrium states. 

One could argue that numerical models often at best conserve energy and enstrophy, 
so there would be at first sight no reason to consider higher order invariants to predict 
the simulated self-organization. However, the truncated dynamics in numerical models 
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may lead to the formation of quasi-stationary states having a life time that tends to 
infinity as the resolution is refined. These quasi-stationary states are distinct from the 
energy-enstrophy equilibria, and thus require more invariants to describe using a statisti- 
cal theory. The formation of such quasi- s tation ary states is a characteristic of long-range 
interacting systems (e.g., ICampa et all (2009) and references therein). Given these re- 
sults, our second working hypothesis is that the quasi-stationary states found in the 
truncated dynamics of numerical simulations may be close to a state predicted by MRS 
theory. 

There are only a few studies of s t atistic al equilibria of the continuously stratified quasi- 
geostrophic dynamics. iMerrvfieldl l| 19981) obtained and discussed a class of solutions in 
the framework of energy-enstrophy statistical theory above topography, for which there is 
linear relation between potential vorticity q and streamfunction tp. The solutions obtained 
were critical points of the theory, but the actual entropy maxima were not selected. Here 
we consider the particular class of MRS equilibrium states characterized by a linear q — ip 
relation, so we describe the same class of equilibrium states, except that we do not take 
into account bottom topography. The effect of bottom topography will be addressed in 
a future work. The novelty of our work is to relate each critical point to the energy and 
fine-grained enstrophy constraints and to select the actual entropy maxima in different 
limit cases. 



Schecterl (|2003l ) provided a direct comparison between numerical simulations and nu- 



merical computation of MRS statistical equilibria, focusing on the particular case of 
initial conditions given by a random field statistically homogeneous in space, as in 
iMcWilliams et all (|1994h . This corresponds to a particular case for which the enstro- 
phy profile is only weakly varying on the vertical. We will outline in this paper the 
important consequence of the conservation of enstrophy at each depth. This important 
role played by the enstrophy constraint allows us to discuss how the beta effect can 
favor barotropization. For a given initial streamfunction field, the beta effect modifies 
the initial distribution of potential vorticity, and therefore modifies the vertical profile 
of fine-grained enstrophy. More precisely, increasing the beta effect tends to increase the 
depth independent part of the fine-grained enstrophy profile. We will show that sta- 
tistical equilibria associated with depth independent fine-grained enstrophy profiles are 
barotropic, which means that increasing the beta favors barotropization, according to 
statistical mechanics predictions. 

The paper is organized as follows: the continuously stratified quasi-geostrophic equa- 
tions are introduced on section [2j The MRS equilibrium statistical mechanics theory of 
such systems is presented in section 3. Here it is shown that the class of MRS statistical 
equilibria characterized by a linear relation between potential vorticity and streamfunc- 
tion are solutions of a minimum enstrophy problem, which is solved in various limit cases, 
first without beta effect and then with beta effect. The theory predicts an increase of 
barotropization with increasing value of beta, which is tested against numerical simula- 
tions in section 4. A summary of the main results and a discussion of their application 
to the ocean is given section 5. 



2. Continuously stratified quasi-geostrophic model 

2.1. Equations of motion 

We consider an initial value problem for a freely evolving, unf orced, undissipated, quasi- 
geostrophic flow with continuous stratification N(z) (see e.g.. IVallia (|2006l ). section 5.4): 
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d t q + J(i>,q)=0, (2.1) 



with J(tjj,q) = d x tpd y q — d y ij)d x q , and / = / + j3y is the Coriolis parameter with the 
/3-plane approximation. Neglecting buoyancy and topography variations, the boundary 
condition at the bottom (z = —H, where H is a constant) is given by 

d z ^\ z =- H = 0. (2.3) 

The boundary condition at the surface (defined as z — 0, using the rigid lid approxima- 
tion) , is given by the advection of buoyancy 

d t b s + J($,b s ) = 0, ^0,V|*=O = &.- (2-4) 

The reason we focus on quasi- geostrophic dynamics is that f or these f lows, there is 
no direct cascade of energy and no dissipative anomaly, see e.g. I Salmon (1998), which 



makes possible a straightforward generalization of statistical mechanics theories devel- 
oped initially in the framework of two-dimensional Euler equations. As far as the ocean 
is concerned, quasi-geostrophic dynamics is relevant to describe "mesoscale turbulence", 
but does not account for "sub-mesoscale turbulence", which in volve ageostrophic ef - 



fects, with part of the en ergy that may cascade forward, see e.g. iThomas et al\ (2008); 
Ferrari fc Wunsch (2009) and references therein. 



2.2. Modal decomposition 

2.2.1. Laplacian and baroclinic eigenmodes 

We consider the case of a doubly periodic domain T> = [— n n] x [— 7r n]. The Fourier 
modes ek,i(x,y) are a complete orthonormal basis of the 2D Laplacian Aefc,; = —K 2 ek,i 
where (fc, I) are the wavenumbers in each directions and K 2 = k 2 + l 2 . In the vertical, the 
complete orthonormal basis of barotropic-baroclinic modes is defined by the solutions 
F m (z) of the Sturm-Liouville eigenvalue problem 

Tz {^z Fm ) = " A - Fm ^ ( 2 - 5 ) 

where m > is an integer, Aj^ is the eigenvalue associated with mode m, which defines 
the m'^-baroclinic Rossby radius of deformation A^ 1 . The barotropic mode is the (depth 
independent) mode F , associated with A = 0. 

2.2.2. Surface quasi-geostrophic (SQG) modes 

In the following, it will be convenient to distinguish the interior dynamics from the 
surface dynamics. For a given potential vorticity field q(x,y, z,t) and surface buoyancy 
field b s (x,y,t) , the streamfunction can be written ip = ipi nt + Vwf, where ipi nt is the 
solution of 

^ + 4- (tStttV) =q-Py with 0,^U=o = 3^U=_h = o, (2.6) 



dz \N 2 dz 
and where VWf is the solution of 
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Each projection of i/wf on a (horizontal) Laplacian eigenmode (I, k) defines a surface 
quasi-geostrophic eigenmode (SQG mode hereafter). When the stratification N is con- 
stant and when fo/(NK) H, these SQG modes are decreasing exponential functions 



(2.8) 



where the coefficients Ak,i are determined using the boundary condition at the surface. 
When the stratification is not constant, but when fo/N(0)K <C H, the SQG modes are 
still characterized by a decreasing exponential, with an e-folding depth hj< = fo/KN(0). 

The distinction between surfa ce and interior streamfu nction is sometimes useful to 
describe the dynamics , see e.g. Lapevre &: Klein (|2006l) . but a shortcoming of such a 



decomposition is that ^;„t and Vwf are not orthogonal. 

2.3. Delta-sheet approximation and SQG-like modes 

The boundary condition (|2.4|) can be formally replaced by the condition of no buoyancy 
variation (d z ip = at z = 0), provided that surface buoyanc y anomalies are in terpreted 
as a thin sheet of potential vorticity just below the rigid lid ([Brethertonl 1 1 9661 ) . For this 
reason, and without loss of generality, we will consider that b s = in the remainder of 
this paper. 

Let us now consider a case with no surface buoyancy, but with a surface intensified 
potential vorticity field defined as q = q SU ii(x, y)Q{z + Hi) + qi n t{x, y, z), with Hiq sur { 
-fffftnt, H% <C H and O the Heaviside function. Using the linearity of (|2.6[) . one can still 
decompose the streamfunction field into ip — Vwf + V^nt where V'surf is the flow induced 
by the PV field q — q SUI {Q(z + Hi), obtained by inverting (|2.2p . and tpi nt the flow induced 
by 9int- For z < —Hi, the streamfunction Vwf satisfies equation (|2.7|) . just as SQG 
modes. This is why we call it an SQG-like mode. 



3. Equilibrium statistical mechanics of stratified quasi-geostrophic 
flow 

In this section we briefly introduce the MRS equilibrium statistical mechanics theory, 
and introduce a general method to compute the MRS equilibria. This allows for explicit 
computation of statistical equilibria on the /—plane in subsection [331 an d on the /3-plane 
in subsection 13.41 

3.1. Theory 

Let us consider an initial (fine-grained) PV field qo(x,y,z). The dynamics (|2.1j) induced 
by this PV field stirs the field itself, which rapidly leads to its filamentation at smaller 
and smaller scale, with the concomitant effect of stretching and folding. There are two 
(related) consequences, one in phase space, one in physical space. 

In physical space, the filamentation toward small scales implies that after a sufficiently 
long time evolution, the PV field can be described at a macroscopic level by a probability 
field p(x, y, z, a, t)da to measure a given PV level between a and a + da at a given point 
(x, y, z). This probability distribution function (pdf hereafter) p is normalized locally 

Mx,y,zN[p]= I da p = 1, (3-1) 

where the integral is performed over the range S of PV levels, which is prescribed by 
the initial condition. The field p is a macrostate because it is associated with a huge 
number of fine-grained PV configurations q(x, y, z) (the microscopic states), in the sense 
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that many realizations of the (microscopic, or fine-grained) PV field q leads to the same 
macroscopic state p. 

The phase space is defined by all configurations of the field q(x, y, z). The cornerstone 
of the theory is to assume that stirring allows the system to explore evenly each possible 
configuration of the phase space that satisfies the constraints of the dynamics, so that 
these configurations can be considered equiprobable. When this assumption fails, as in 
the linear regime discussed in section 4.4, the statistical theory fails. 

The main idea of the theory is then to count the number of microstate configurations 
associated with a macrostate (the pdf field p). and to say that the equilibrium state is 
the most probable macrostate that satisfies the constraints of the dynamics, i.e. the one 
that is associated with the largest number of microstates satisfying the constraints. One 
can then show that a n overwhelming number of microstates are associated with this most 
probable macrostate ( Michel fe Robertlll994 ). An important physical consequence is that 



an observer who wants to see the actual equilibrium macrostate would neither need to 
perform ensemble average nor time average, but would simply have to wait a sufficiently 
long time. 

Note that the output of the theory is the pdf field p, but the quantity of interest to 
describe the flow structure is the coarse-grained PV field defined as 

q= da a p. (3.2) 

In practice, counting the number of microstates q associated with a macrostate p is a 
difficult task, b ecause q is a continuo u s field , and it requires the use of large deviation 
theory, see e.g. Bouchet fc Corvellec ( 2010h . However, it can be shown that classical 



counting arguments do apply here, and that the most probable state p maximizes the 
usual Boltzmann mixing entropy 

S = — dxdy / dz I da p\np , (3-3) 



V J-H JT. 



sec 



Bouchet fc Corvelled (|2010h and references therein for further details. 



Let us now collect the constraints that p must satisfy, in addition to the local normal- 
ization (|3.ip . An infinite number of constraints are given by conservation of the global 
distribution of PV levels Via, z), which is prescribed by the initial global distribution of 
PV levels P (cr,z): 



V(a,z)[p} = / dxdy p = P {a,z). (3.4) 
Jv 

These conservation laws include the conservation of the fine-grained enstrophy 



dxdy / da a" p = Zq . (3-5) 
Another constraint is given by the conservation of the total energy 

Because an overwhelming number of microstates are close to the equilibrium state p, the 
energy of the fluct uations are negligible w ith respect to the energy due to the coarse- 



grained PV field q ([Michel fc Robertlll994l ). Integrating [3T6] by parts and replacing q by q 
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gives the constraint Eq — £[p], where Eq is the energy of the initial condition, and where 
£=- I dz I dxdy (f3y - q) = - j dz j dxdy j da {fly - a) tpp . (3.7) 

Finally, the MRS theory provides a variational problem 

S(Eo,Po(<r,z))= max {S [p] \ £ [p] = E Q k V [p] = P (<t, z)} , (3.8) 

which means that the equilibrium state is the density probability field p Tsm that satisfies 
the local normalization constraint (|3.1I) . the energy constraint £ [p] = Eq with £ given 
by ()3.7j) . the incompressibility constraint V{a 1 z) [p] — Po(a, z) with V given by p.4j) . 
and that maximizes the entropy functional S given by (|3 . 3[) . We see that the parameters 
of this problem are the values of the constraints, which are prescribed by the initial 
condition q§{x, y, z). 

3.2. MRS equilibrium states characterized by a linear q — ip relation 

3.2.1. Simplification of the variational problem 

In practice, MRS equilibrium states are difficult to compute, because such compu- 
tations involve the resolution of a variational problem with an infinite number of con- 
straints. Computing the critical points of (|3.8I) is straightforward, see e.g. equation (|A 1[) 
of appendix A, but showing that second order variations of the entropy functional are 
negative is in general difficult. However, recent results have allowed considerable s impli- 
fications of the analytical computations of these equilibrium states ( Bouchetll2008 ). The 



idea is first to prove equivalence (usually for a restricted range of parameters) between 
the complicated variational problem f|3 . 8[) and other variational problems more simple 
to handle, involving less constraints than the complicated problem, but characterized by 
the same critical points, and then to explicitly compute solutions of the simpler varia- 
tional problem. The equivalence between the variational problems ensures that solution 
of the simple problem are also solutions of the more complicated one. We will follow this 
approach in the following. 

3.2.2. Minimum enstrophy variational problem 

Starting from the MRS variational (|3.8p . it is shown in appendix A that any MRS 
equilibrium state characterized by a linear q — ip relation is a solution of the variational 
problem 

Zlf {E , Zo(z)) - mm {Zg [q] \£[q] = E } (3.9) 

with 

Zlf [q}= [ dz^L, Z cg = \\ dxdy q\ (3.10) 

J-H Z ^ JT> 

where Zq{z) is the fine-grained enstrophy profile, Eq the energy. This variational problem 
amounts to find the minimizer q m (a coarse grained PV field) of the coarse-grained 
enstrophy (with [q m ] = {Eq, Zq{z))), among all the fields q satisfying the 
constraint £\q] = Eq given by (|3.6j) . 

Critical points of the variational problem (|3.9p are computed by introducing a Lagrange 
multiplier f3 t associated with the energy constraint, and by solving 

szlf + we = o, 

where first variations of the functionals are taken with respect to the field q. Because 
fit is the Lagrange parameter associated with the energy constraint, it is called inverse 
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temperatur^j]. Using SZ^* = J v dxdy f® H dz qSq/Zo and S£ — — J dxdydz ipSq, one 



finds that critical points of this problem are 

q = PtZoi). 

The flow structure of these critical points can then be computed by solving 



(3.11) 



A ^+7TTkir^ + f 3 y = Z °i 3 ^ with dzMz=o,*=-H = 0. (3.12) 
az N z az 

Critical points can then be classified according to the value of their inverse temperature 
fit ■ All the problem is then to find which of these critical point are actual minima of the 
coarse-grained enstrophy [q] for a given energy £ \q] = Eq . 

MRS equilibrium states characterized by lin ear q — ij) relations are expe cted either in a 
low energy limit or in a strong mixing limit, see Bouchet fc Venaille] ( 2012 ) and references 
therein. Let us estimate when the low energy limit is valid. For a given distribution 
Po(a,z) of PV levels a, one can compute the maximum energy E max among all the 
energies of the PV fields characterized by the distribution Pq{<j, z). Low energy states 
characterized by the same distribution Pq{<t, z) are therefore those characterized by E <C 
Emax- The strong mixing limit, which can in some case overlap the low energy limit, 
is obtained by assuming f3 t aip <C 1 in the expression (IA II) of t he critical points of the 
MRS variational problem (|3.8[) . see ( Chavanis fc Sommerialll996[ ). Note that these results 
would remain valid in the presence of bottom topography. The only change would be the 
expression of the energy, and the bottom boundary condition. 

3.2.3. Link with the phenomenological principle of enstrophy minimization and other 
equilibrium states 

At a phenomenological level, the variational problem (13.911 c an be interpreted a s a 
generalization of the minimum enstrophy principle of iBretherton fc Haidvogel ( 19761 ) to 
the continuous stratified flow: the equilibrium state minimizes the vertical integral of the 
coarse-grained enstrophy normalized by the fine-grained enstrophy in each layer, with a 
global constraint given by energy conservation. This normalization is necessary to give a 
measure of the degree o f mixing comparable from one d epth z to another. 

It has been shown by Carnevale fc Frederiksenl dl987l) that solutions of the minimum 
enstrophy problem of IBretherton fc Haidvogell (|1976l ) for one layer quasi-geostrophic 
flows, are energy-enstrophy equilibrium states of the Galerkin-truncated dynamics of 
these models, in the limit when the high-wave number cut-off goes to infinity. These 
energy-enstrophy equ ilibrium states we r e first introdu c ed in the framework of 2D trun- 
cated Euler flows by iKraichnan ( 1967 ) . Salmon et al. dl976l) generalized th e theor y to 
one or few layers truncated quasi-geostrophic flows. iFrederiksen fc Sawfordl (|l980l) ap- 
plied the theory for a barotropic flow on a sphere, and IFrederiksen! (|19916l ) generalized 
these results for two-layers quasi-geostro phic flow on a sphe re. In this latter case, inter- 
esting stability results were obtained by IFrederiksen! ( 1991a ) . 

The case of continuously str atified flow above topography was addressed by Merrvfieldl 
( 1998) for truncated dynamics. Merrvfie idl (|l998l) computed critical points of the statisti- 
cal mechanics variational problem, found q = fi t (z)ip and discussed the vertical structures 
of these flows depending on fi t . We showed in the previous subsection that the parameter 
/i t is related to the fine-grained enstrophy profile through fj, t — f5tZ§(z). In addition, the 



f The notation /3 is traditionally used in thermodynamics for inverse temperature. We use 
here the convention fit for this inverse temperature, since j3 = d y f refers to the beta effect, i.e. 
to the variation of the Coriolis parameter with latitude. 
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critical states described by iMerrvfieldl (jl998h could either be saddle, minima or max- 
ima of the entropy functional. We show in the next section that we can select which of 
them are actual entropy maxima (or equivalently coarse-grained enstrophy minima), in 
different limit cases. 

3.2.4. Transients, mean and sharp equilibrium states 

In the framework of energy-ens trophy theory, a disti n ction is often made between "tran- 
sients" and "mean" energy, see e.g. lSalmon et al. l (|l976h ; lFrederiksed(fl99Tl ; lMaida fc Wand 
(120061) . among others. Statistical equilibria are called "sharp" when the "transients" van- 
ish. In earlier works, the "mean" quantities were often computed by summing over all 
the microstates distributed according to the canonical measure, which were implicitly 
assumed to be equivalent to the microcanonical measure. In the absence of topography 
and beta plane, the "mean" states defined in this way are always zero in a doubly peri- 
odic domain, whatever the resolution of the truncated model, and all the energy is in the 
"transients". 

By contrast, the equilibrium states computed in the MRS framework in the micro- 
canonical ensemble are always sharp: the energy of t he fluctuations necessarily va nish, 
which can be shown with large deviation theory, see iBouchet fc Corvelled (|201dh and 
references therein. This apparent paradox is solved when one realizes that microcanon- 
ical and canonica l ensembles are, i n general, not equivalent for lon g-ranged interacting 
system s, see e.g. Ellis et al. (2000) for general considerations and Venaille fc Bouchet 
(|2011a ) for application to two-dimensional and geophysical flows. 

In the case of a doubly periodic domain without bottom topography and beta effect, 
the MRS equilibrium states (or the e nergy-enstrophy equilibria) are degenerate, because 
of the symmetries of the domain, see (|Bouchet fc Venailldl2012li for more details and for 
a discussion on how to remove the degeneracy, see also lHerbert et al. ( 20121 ) for the case 
of a spherical domain. In practice, the dynamics breaks the system symmetry by selecting 
one of the degenerate equilibrium states. Of course, if one considers the average potential 
vorticity field over all the degenerate equilibria, then one recovers zero, and this is what 
happens when computing the "mean" state in the canonical ensemble. We stress that for 
either quasi-geostrophic flows or truncated quasi-geostrophic flows in the limit of infinite 
resolution, the freely evolving dynamics itself can not jump from one state to another 
once the symmetry is broken: in the microcanonical ensemble, there are no "transients". 
We also stress that the physically relevant statistical ensemble for an isolated system 
such as a freely evolving flow is the microcanonical ensemble. 

3.3. Enstrophy minima on a f -plane 
In the previous subsection, we have shown that MRS equilibrium states characterized 
by a linear q — ip relationship are solution of a minimum coarse-grained enstrophy varia- 
tional problem (|3.9p . Such MRS equilibrium states will be referred to as "coarse-grained 
enstrophy minima" in the following. We found in the previous section that the affinc 
relation of these coarse-grained enstrophy minima are of the form q = fStZgip, with a 
coefficient proportional to the fine-grained enstrophy Zq, and are solutions of (|3.12[) . The 
critical points of the variational problem (|3.9j) are any states q satisfying these condi- 
tions. The aim of this subsection is to find which of these critical point q are the actual 
coarse-grained enstrophy minima, solutions of variational problem (|3.9p . when there is 
no f3— effect, but for arbitrary stratification N. 

Injecting ()3.11[) in (|3.10ll . and using the expression ()3.7|) with the constraint Eq — £[p], 
one finds that the coarse-grained enstrophy of each critical point q is given by 



-PtEo ■ 



(3.13) 
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N 




Figure 1. Left panel: three different fine grained enstrophy profiles. Right panel: corresponding 
vertical structure of statistical equilibrium states (if>(z)/ip(0) on the left panel), in the case of 
constant stratification ( fo/N 2 — 0.1). The e-folding depth in case A is h = fo/NK (with here 
K = 1 for the statistical equilibrium state). 

Remarkably, the coarse grained enstrophy of a critical point q depends only on the 
inverse temperature /3 t . We conclude that coarse-grained enstrophy minimum states are 
the solutions of (|3.12[) associated with the largest value (3f 
Projecting (|3.12p on the Laplacian eigenmode ei t k{x,y) gives 

!(jy^w) = (M + J^w. (3-14) 

with 

d z 4>k,i\z=o,-H = 0, K 2 = k 2 + l 2 , tp = ^2$k,ieki- 

k,l 

We see that each critical point is characterized by a given wavenumber modulus K. Its 
vertical structure and the corresponding value of /3 t must be computed numerically in 
the case of arbitrary profiles Zo(z). Let us consider the example shown in Fig. [TJ for a 
two-step fine-grained enstrophy profile 

Z = Z maf Q(z + H 1 ) + Z bjt Q(-z-Hi), H X 4HH, (3.15) 

where O is the Heaviside function, and for Z- m t/Z svlI { varying between and 1. We find 
that the minimum coarse-grained enstrophy states are always characterized by the gravest 
horizontal mode on the horizontal (K — 1). As for the vertical structure, we observe Fig. 
Q]a tendency toward more barotropic flows when the ratio Zi nt /Z snT { tends to one. 

It is instructive to consider the limit cases Z SUT { = Z- mt and Z- mt — 0, corresponding 
respectively to case A and C of Fig. [TJ for which minimum enstrophy can be explicitly 
solved. When the enstrophy profile is depth independent (Z (z) = Z int , case C of figure 
[TJ, solutions of (13. 14)) are given by the modes F m (z)eK(x,y) defined in subsection (I2.2[) . 
and are associated with inverse temperatures f3 t = — (A^ + K 2 } /Z- mt . We see that the 
maximum value of (3t is reached for the gravest horizontal mode (K = 1) and the gravest 
vertical mode (Ao = 0). It means that the coarse-grained enstrophy minimum state is 
barotropic. 

When Zi ut = (case A of Fig. [TJ, it is straightforward to show that critical points, 
the solutions of ()3.14j) . are SQG-like modes (see subsection 12. 3j) of e-folding depth h = 
fo/N(0)K, associated with fit — —N(0)K/(HiZ SUT ifo) with corrections of order H\/H 
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and of order h/H. The coarse-grained enstrophy minimum is therefore the SQG-like 
mode associated with K — 1 and with the largest e— folding depth h — Lf /2irN, where 
L is the domain length scale. 

These examples show the importance of the conservation of fine-grained enstrophy to 
the vertical structure of the equilibrium state. The main result is that statistical mechan- 
ics predicts a tendency for the flow to reach the gravest Laplacian mode on the horizontal. 
The vertical structure associated with this state is fully prescribed by solving Q3.14p with 
K = 1. Because the barotropic component of such flows are larger than solutions of (|3.14l) 
with K > 1, we can say that the inverse cascade on the horizontal is associated with 
a tendency to reach the gravest vertical mode compatible with the vertical fine-grained 
enstrophy profile Zq. This means a tendency toward barotropization, although in general, 
the fact that the profile Zq is non constant prevents complete barotropization. 

Note that the previous results are presented in the case of a doubly periodic domain, 
but it would be straightforward to generalize them to channel geometries, or to any 
other domains with boundaries. The only change would be the spatial structure of the 
Laplacian eigenmodes. For these domain geometries, one would also need to discuss the 
effect of the circulation values at the boundaries. 

3.4. Including the j3- effect . 

We now discuss how the results from the previous section generalize in the presence of 
beta. For a given initial condition ipo(x,y, z), increasing fj increases the contribution of 
the (depth independent) available potential enstrophy defined as 



to the total fine-grained enstrophy profile Zq(z) = J v dxdy q$, where qo is the initial PV 
field that can be computed by injecting %j)Q in (|2.2j) . For sufficiently large values of f3, the 
PV field is dominated by the beta effect (qo » f3y), Zq therefore tends to Z p and becomes 
depth independent. Because statistical equilibria computed in the previous subsection 
were fully barotropic when the fine-grained enstrophy Zq was depth- independent, we 
expect a tendency toward barotropization by increasing fj. 

Consider for instance the case where tpo is an SQG-like mode such that go — fiy f° r 
z < —Hi, with Hi <C H and associated with a surface fine-grained enstrophy Zq = Z SUT i 
for z > —Hi. The initial profile of fine-grained enstrophy Zq is thus given by (|3.15|) . with 
the interior enstrophy Z lrA — Z p prescribed by the available potential enstrophy (|3.16[) . 
The question is then to determine whether the results obtained previously without beta 
effect still hold in presence of it. 

In oder to apply properly MRS theory in the presence of beta, one must consider 
southern and northern boundaries. Indeed, MRS theory relies on the conservation of the 
global fine-grained PV distribution. Strictly speaking, this distribution is not conserved 
in the doubly periodic configuration on a beta plane, since any fluid particle that travels 
2n in the y direction gains a value — 2tt/3, and one can not know a priori how many 
times a fluid particle will turn over the domain in the y direction before the flow reaches 
equilibrium. In practice, the theory can still be used in the doubly periodic case to make 
qualitative predictions, as explained in the next section. 

In the channel configuration, the streamfunction is a constant at y — ±tt, and we 
can compute the statistical equilibria in two limit cases. First, when fj — 0, we recover 
the configuration with surface intensification of the fine grained enstrophy profile Zq, 
for which we found that statistical equilibria were the SQG-like modes associated with 
the gravest Laplacian horizontal modes, see subsection 13.31 Second, with fj sufficiently 




(3.16) 
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large, such that Z p 3> Z Bur { , the statistical equilibria can be computed by considering the 
case of a constant fine-grained enstrophy profile Zq. In the general case, the values of the 
streamfunction at the northern and southern boundaries are different, and are determined 
using constraints given by mass conservation and circulation conservation along each 
boundary, see e.g. (|Pedloskvl[l982h . Here we consider the simple case ipix, y = ±7r, z) = 0, 



with zero total circulation (J v dxdyq = 0) at each depth. It is necessary to take into 
account the conservation of linear momentum (associated with translational invariance 
of the problem in the x direction), which provides an additional constraint C\q] = Lq, 
with 



m = f < 

J-H 



dxdy qy. (3-17) 



v 



Critical points of the variational problem (|3.9p with this additional constraint satisfy 

q=(M-HV)Z , (3.18) 

where /i and fit are respectively the Lagrange multipliers associated with the linear 
momentum conservation and the energy conservation. Equilibrium states therefore satisfy 

Z (M - W ) = + ^ (|jf ^V>) + ■ (3-19) 

Taking /i = —fi/Zo, we recover the /-plane case of the previous subsection: solutions of 
Eq. (|3.19p are barotropic-baroclinic modes on the vertical with eigenvalue — A^, and are 
Laplacian eigenmodes in the channel (with eigenvalue denoted by —K 2 ), with /3 t Z = 
—K 2 + X 2 n . These states characterized by fi t = — (K 2 + A^) /Z , /i = —f3/Z are drifting 
towards the west at speed of Rossby waves: Udrift = -PI (>4n + R2 )- The enstrophy 
minimizer among these states is the one associated with the smallest value of fit > namely 
the barotropic mode associated with the gravest Laplacian eigenmode consistent with the 
circulation and momentum constraints. When /i ^ P/Zq, Eq. (|3.19|l can be inverted, and 
ip is barotropic since (fi — /iZo)y is depth-independent. The computation of the minimum 
enstrophy state is therefore equivalent to computing the minimum enstrophy state of 
a barotropic flow in a channel. A co mplete discussion of mini mum enstrophy s t ates i n 



a barotropic channel are presented in Corvellec ( 2012 ). see also Maida &: Wane (2006), 



with detailed computations including the effect of bottom topography and of non-zero 
circulations at the boundaries. The important point here is that in the large beta limit, 
the statistical equilibrium states are barotropic. 



4. Numerical simulations 

4.1. Experimental set-up 

We consider in this section the final state organization of an initial SQG-like mode 
(defined in subsection I2.3[) . varying the values of fi in a doubly periodic square domain 
of size L = 2tt. More precisely, we consider a vertical discretization and an horizontal 
Galerkin truncation of the dynamics (|2.1j) . for an initial potential vorticity field go = 
1o,surf(x, y)9 (z + Hi) + jiy, such that go = Py in the interior (— H < z < —Hi) and 
qo ~ qosurf in a surface layer z > —Hi. The surface PV q S urf(x,y) is a random field with 
random phases in spectral space, and a Gaussian power spectrum peaked at wavenumber 
Kq = 5, with variance SKq = 2, and normalized such that the total energy is equal to 
one (Eq = 1). This corresponds to the case discussed subsection (13. 4p . with a vertical 
profile of fine-grained enstrophy given by (|3.15|) , where the interior enstrophy is given by 
the available potential enstrophy (|3.16|l : for z < —Hi, Zq = 47r 4 /3 2 /3. 
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We perform simulations of the dynamics by considering a vertical discretization with 10 
layers of equal depth, horizontal discretization of 512 2 , H = 1, and F = (Lfo/HN) 2 = 1, 
using a pseudo-spectral quasi-geostrophic mo del for the horizontal dynamics without 



small scale dissipation ( Smith fc Vallisl feoOlF ). We choose Hi = H/10 for the initial 



condition, so that there is non zero enstrophy only in the upper layer in the absence 
of a beta effect. We also perform experiments in the presence of small scale dissipation 
(Laplacian viscosity or hyperviscosity), and found no qualitative differences as far as the 
large scale horizontal structure and the vertical structure of the flow were concerned. 
The only difference in presence of small scale dissipation is the decay of the enstrophy 
in each layer, because small scale dissipation smooths out PV filaments. Because of the 
energy cascade, the energy remains constant at lowest order even in presence of small 
scale dissipation. 

Time integration of the freely evolving dynamics proceeds for about 250 eddy turnover 
times. Typically, i) the unstable initial condition leads to a strong turbulent stirring and 
the flow self-organizes rapidly into a few vortices, which takes a few eddy turnover times 
ii) same sign vortices eventually merge together on a longer time scale (a few dozen 
of eddy turnover times), and iii) the remaining dipole evolves very slowly on a slow 
time scale (more than hundreds of eddy turnover times). In absence of /3-effect, a good 
indicator of the convergence is given by the convergence of the q — ip relation to a single 
functional relation. As we shall see in the following, this phenomenology is complicated 
by the presence of beta. 

We saw in the previous section that a proper application of the MRS theory in presence 
of a beta plane would require considering northern or southern boundaries, since the fine- 
grained PV distribution is not conserved in the doubly-periodic domain case. However, 
we still can use the theory to interpret qualitatively the simulations in the doubly periodic 
case. For small values of (3, planetary vorticity gradients are negligible with respect to 
relative vorticity and stretching terms in the PV expression, in which case one can check 
that the fine grained PV distribution remains nearly the same during the flow evolution, 
and MRS can be applied. For large values of (3, the fluid particles can not escape in the 
meridional direction, they are confined in zonal bands having a width of the order of the 
Rhines scakj}] Lr = 2ir ~ E^ 4 / /3 1 / 2 , and thus the fine-grained PV distribution remains 
close to the initial one, and MRS theory can be applied. 

4.2. SQG-like dynamics (f3 = 0) 

The case with f3 = is presented in Fig. [5] The initial PV and streamfunction fields have 
horizontal structures around wave number Kq ~ 5, which is associated with an SQG-like 
mode of typical e-folding depth fo/(NKo) w 0.1. 

An inverse cascade in the horizontal leads to flow structures with an horizontal wavenum- 
ber K decreasing with time, associated with a tendency toward barotropization: the 
e-folding depth of the SQG-like mode increases as fo/2NK. According to statistical me- 
chanics arguments of the previous section, the flow should tend to the "gravest" horizontal 
SQG-like mode (K = 1), which is also the "gravest" vertical SQG mode, i.e. the one as- 
sociated with the largest e-folding depth. The concomitant horizontal inverse cascade 
(most of the kinetic energy is in the gravest horizontal mode K = 1 at the end of the 
simulation) and the increase of the e-folding depth are observed on Fig. [5J showing good 
qualitative agreement between statistical mechanics and numerical simulations. 

There are however quantitative differences between the observed final state organiza- 

1/2 

f Here Eo = 1 is the initial energy, so that E is a good metric for typical velocities of the 
flow (assuming it becomes barotropic). 
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tion and the prediction of MRS equilibria characterized by a linear q — ip relation. As 
seen on Fig.[3J the q — ip relation is well defined, but it is not linear. Taking into account 
higher order inv ariants in the MRS theory w ould allow to better describe this sinh-like 
relation, see e.g. ( Bouchet fc Simonnet 2009), but it may not be sufficient to describe the 



"unmixed" cores of the remaining vortex in the final dipole (jSchect er 2003) . 

The presence of these "unmixed" PV blobs implies the existence of structures smaller 
than wavenumber K = I. Because each projection of the PV field on a different wavenum- 
ber is associated with an SQG-like mode associated with a different e-folding depth, the 
vertical structure of the flow is actually a linear combination of SQG-like modes associ- 
ated with different wavenumbers. A consequence is that the effective e-folding depth of 
the kinetic energy is fo/(2NK cS ). We estimated the coefficient K e g as 3.5 by considering 
the linear q — ip relation passing through the extremal points of the observed q — ip relation 
of Fig. [3] Let us conclude on this case by discussing the relaxation towards equilibrium: 
we observed that the q — ip relation did not change at all between 150 and 250 eddy 
turnover times, proving that the system reached a stationary state. In addition, we see 
in Fig. [SJ-a that the kinetic energy of each vertical mode (and therefore of each layer) 
reaches a plateau after a few eddy turnover times. As expected from the MRS theory, 
the statistical equilibrium is sharp, there is no transient. 



4.3. Effect of beta (/3 ~ 1 ) 

We now consider the free evolution of the same initial condition ipo(x, y, z) as in the pre- 
vious subsection, but on a beta plane. Because ft ^ 0, the initial condition for potential 
vorticity is different than in the previous section: the total energy and initial velocity 
fields are the same as before, but the vertical profile of fine grained enstrophy is dif- 
ferent, since there is now available potential enstrophy in the interior. As explained in 
subsection 13. 3[ this contribution from interior fine-grained enstrophy shall increase the 
relative barotropic component of the statistical equilibrium state. This is what is actually 
observed in the final state organization of Fig. [4] We conclude that in this regime, the 
beta effect is a catalyst of barotropization, as predicted by statistical mechanics. 

Additionally, the q — ip relation presented in Fig. [5] presents a very good agreement 
with the prediction of our computation of barotropic MRS equilibria carried in subsection 
13.41 we observe a linear relation between (q — j3y) and ip with a slope faZo = — 1. This 
means that the dipole structure observed in the streamfunction field is the actual gravest 
horizontal Laplacian eigenmode K = 1, and that this eigenmode is drifting westward at a 
speed (3/K 2 . In addition, we see in Fig.[8]that after a few eddy turnover times, almost all 
the kinetic energy is in the barotropic mode, and that the energy levels in each vertical 
mode reach a plateau: there are no transients. 

We conjecture that the good agreement with MRS theory occurs since the initial 
distribution of PV levels and the initial energy place the flow in the low energy regime 
considered in the Appendix. In addition, for intermediate values of /3, the Rhines scale 
Lr = 2ttEq 4 I (3 1 / 2 is of the order of the domain scale. The beta effect limit the meridional 
displacements of fluid particles at this scale; this is why the computations performed 
Appendix 2 in a channel geometry are qualitatively relevant to this case. 

4.4. Effect of strong Rossby waves (/3 ^> 1) 

For large values of beta there is still an early stage, lasting a few eddy turnover times, 
of turbulent stirring of the unstable initial condition. But this initial stirring is limited 



t By catalyst we do not mean a dynamical effect, as in chemistry, but simply the fact that 
barotropization in the final state is enhanced in presence of a beta plane. 



The catalytic role of beta effect in barotropization processes 15 



Potential vorticity q g Streamfunction ¥ KE(z) 




Figure 3. q — ip relation in the upper layer for t — 250 eddy turnover time. 
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Figure 4. case f3 = 2.1. See Fig. [2]for legend. Note that the colors scales are changed at each 
time; this is why the interior " beta plane" is not visible in the upper panel of potential vorticity. 



to scales smaller than the Rhines scale Lr = 2t7E^ 4 / ft 1 / 2 , which eventually leads to a 
zonal PV field perturbed by Rossby waves around this scale, see Fig. [6l The dynamic s is 
then dominated by Rossby waves, just as in a single layer model iRhinea (|l977l . ll97.^ . 

Because the perturbation is confined in the upper layer, irreversible turbulent stirring 
occurs in the upper layer only, but not at the bottom. We did check that the global 
distribution of coarse-grained PV levels did not change in the lower layer through time, 
while this distribution was changed in the upper layer. In this case, the fundamental 
assumption that the system explores evenly the available phase space through turbulent 
stirring is not valid. One might reasonably argue that the statistical mechanic theory 
is incomplete, for it cannot identify the phase space that can be explored. In any case, 
we conclude that large values of beta prevent the convergence toward this equilibrium, 
leaving the question of the convergence itself to future work. 
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Figure 5. Relation between q — fly and streamfunction in the upper layer, at t — 100 eddy 
turnover times, for j3 — 2.1. The fact that this relation is close to be linear shows that the 
small energy limit considered in this paper is relevant to this case. The coefficient of the slope is 
approximatively one, meaning that the energy is condensed in the gravest horizontal Laplacian 
eigenmode, which is propagating westward at the speed of barotropic waves. This is consistent 
with the computation of minimum enstrophy states on a beta plane, see subsection 13.41 



4.5. Turbulent stirring vs wave dynamics 

We summarize the effect of planetary vorticity gradients on Fig. [7] by considering the 
ratio E^o/Eq of barotropic energy to total energy, after 250 eddy turnover times, for 
varying values of f3. The plain red line represents statistical mechanics predictions. In 
order to allow comparison with the numerical experiments, for which the horizontal PV 
field did present structures smaller than the gravest mode K = 1 for low values of (3 (see 
subsection l4.2|) . we computed the vertical structure with equation (|3.14p . using K c s = 3.5, 
and using the enstrophy profile Zq(/3) of the initial condition Strictly speaking, this 
effective horizontal wavenumber K c g should vary with beta between K e s = 3.5 for f3 — 
and K e s — 1 for f3 > 1, but taking into account this variations would not change much 
the shape of the red plain curve. 

The critical value of /3 between the turbulent stirring regime for which statistical 
mechanics predictions are useful and the wave regime can be estimated by considering 

the case in which the Rhines scale L^h ~ 2-7T I E Q /p\ is of the order of the domain 
scale L — 2ir. Here the total energy is -Bo = 1> which renders the critical value /3 ~ 1 in 
our simulations. 

Because the ratio E^^m/Eg gives incomplete information on the flow structure, and 
because dynamical information is useful to understand where statistical mechanics fails 
to predict the vertical structure, we represent in Fig. [5] the temporal evolution of the 
energies E m of the baroclinic modetfj] defined in (|2.5[) . 

For small values of beta (see the first panel of Fig. [5]), the dynamics reaches the 
equilibrium after stirring during few eddy turnover times, and the contribution of each 
mode E m simply reflects the SQG-like structure of the final state. For values of beta 
of order one (see the second panel of Fig. [8J , the same stirring mechanism leads to a 
barotropic flow after a few eddy turnover times. For values of f3 larger than one, there 
is still an initial stirring regime, but the duration of this regime tends to diminish with 
larger values of beta. We observe a slow transfer of energy from first baroclinic to the 

f Note that with this notation, E m= o = Eb&ro is the energy of the barotropic mode, which is 
in general different from the total energy of the initial condition Eq. 
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Figure 6. case (3 = 10, see Fig. [2] for legend. 



barotropic mode in the case /3 = 10 (see the third panel of Fig. [5]), but it is unclear 
if it will eventually lead to a condensation of the energy in the barotropic mode. In 
the case (3 = 100 (see the fourth panel of Fig. [8J, there is in average almost no energy 
transfers between the vertical modes, after the initial stirring regime. Again, it is unclear 
if interactions between waves would eventually condense their energy into the barotropic 
mode in a longer numerical run. The strong temporal fluctuations of the contribution of 
each vertical mode visible on the wave regime of figure Fig. [5] may explain why they are 
large variations of / E to t (in addition to the mean tendency of a decay with j3) in 
the wave regime figure [7] 

To conclude, in the turbulent stirring regime (i.e. Lr ~Lor Lr > L) MRS statistical 
mechanics predictions are qualitatively correct, and there are no "transients" (no temporal 
fluctuations) once the flow is self-organized. By contrast, in the wave regime (i.e. Lr <C 
L) , there are strong temporal fluctuations in the mode amplitudes, and we believe 
that statistical mechanics (either MRS of energy-enstrophy equilibria for the truncated 
dynamics) can not account for these observations. 
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Figure 7. Variation of the Barotropic energy Etaro normalized by the total energy Eo, varying 
the ratio of the domain scale with the Rhines scale L/Lr — /3 1 / 2 E ~ 1/ ' 4 'L/2ir, at the end of each 
run, for the same initial condition ipo (here L — 2n, Eo = 1, so only f3 varies). Red continuous 
line: predictions from statistical mechanics. The vertical structure of the equilibrium state is 
computed with equation ()3. 14f) . using K e g = 3.5, and using the enstrophy profile Zo((3) of 
the initial condition ipo. The dashed-dotted line separates the wave regime from the turbulent 
stirring regime, in which we expect the statistical mechanics predictions to be valid. 



5. Discussions and conclusions 

We have used statistical mechanics to predict the flow structure resulting from the 
self-organization of freely evolving inviscid stratified quasi-geostrophic flows. The only 
assumption of the theory is that the flow explores evenly the available phase space. In 
order to compute explicitly statistical equilibria, and to discuss basic physical properties 
of these equilibria, we have focused on MRS equilibria characterized by linear relation 
between potential vorticity and streamfunction. For such MRS equilibrium states, we 
need the total energy and the fine-grained enstrophy in each layer, although all the 
invariants are implicitly taken into account. 

We explained that these states are expe cted in a low e nergy limit or in a strong mixing 
limit. By applying a method proposed by (|Bouchetll2008l) . we have shown that MRS equi- 



libria characterized by a linear q — ip relation are solutions of a minimum coarse-grained 
enstrophy variational problem, which generalizes the Bretherton-Haidvogel minimum en- 
strophy principle to the stratified case. 

The central result of the paper is the elucidation of the physical consequences of the 
conservation of the fine-grained enstrophy Zq(z) on the structure of the equilibrium 
states. We showed first that the affine coefficients between the coarse-grained PV field 
and the streamfunction are proportional to the fine-grained enstrophy: q = (3tZo(z)ip. 
This relation allowed to compute explicitly the statistical equilibria, and to discuss their 
structure depending on the profile Zq(z). When this profile Zo is depth independent, 
the equilibrium state is fully barotropic, whatever the stratification profile N. When this 
profile Zq is surface intensified, the equilibrium state is an SQG-like mode, characterized 
by an e-folding depth ~ fo/NK, where 1/K is bounded by the domain scale. Since 
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Figure 8. Temporal evolution of the contribution of the energy of each baroclinic mode to the 
total energy, for different values of /3. The modes are indexed by m, with m = the barotropic 
mode and m = 1 the first baroclinic mode. 



the e-folding depth increases with K , larger horizontal flow structures imply "more 
barotropic" flows. Since statistical equilibria are associated with the gravest Laplacian 
horizontal modes (K — 1), the equilibrium state is the "most barotropic one" given an 
enstrophy profile Zq. 

We conclude that the dynamics leads to the gravest horizontal mode on the horizon- 
tal (inverse cascade), associated with the gravest vertical structure consistent with the 
conservation of the fine grained enstrophy. The flow becomes fully barotropic when the 
fine-grained enstrophy profile is depth independent. 

These results can be used to understand the role of beta in barotropization in the 
following way: consider an initial SQG-like flow, with /3 = 0, prescribed by its stream- 
function; its enstrophy profile is surface intensified, and so will be the equilibrium sate. 
But when one switches on the beta effect, with the same initial streamfunction, the con- 
tribution of the depth independent part of the fine-grained enstrophy increases, which 
means a tendency toward a more barotropic equilibria, according to the previous conclu- 
sion. This result reflects the fact that in physical space, the initial SQG-like mode stirs 
the interior PV field (initially a beta plane) , which in turn induces an interior flow, which 
stirs even more the interior PV field, and so on. 

When the beta effect becomes large enough, the PV distribution at each depth z 
becomes prescribed by the initial beta plane, so the most probable state should be also 
characterized by a depth independent PV field, i.e. by a barotropic flow. 

To reach the equilibrium state, the PV field must be stirred enough to explore the 
phase space. Yet large values of beta stabilize the initial condition, the system is trapped 
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in a stable state different from the equilibrium state, and the flow dynamics becomes 
dominated by the interaction between Rossby waves rather than by turbulent stirring. 
The reason for this trapping is that the dynamics can not provide by itself sufficient 
perturbations to escape from the stable state and to explore other regions of the phase 
space, which would attract the system toward the equilibrium state. We observed in 
numerical simulations that this wave regime did not lead to a barotropization of the 
flow after a few hundred turnover times. We do not know if still longer numerical runs 
would lead to more barotropization. We estimated the critical value of beta between the 
turbulent stirring regime and the Rossby wave regime as the value such that the Rhines 
scale is of the order of the domain scale, which gives (3 = 4tt 2 E^ 2 / L 2 . 

In this paper, we focused on the effect of j3 on barotropization process. More generally, 
our results show that vertical transfers of energy and momentum are favored by the 
presence of any lateral potential vorticity gradients, since these gradients provide a source 
of available perturbation enstrophy. These lateral PV gradients may for instance be due 
to large scale mean flows set by external process (large scale wind pattern for the oceans, 
and large scale temperature gradients for the atmosphere). A similar situation occurs 
in the framework of interactions between waves (or eddies) and mean flows: horizontal 
momentum can be transfered vertically by inviscid form drag if there are horizontal 
fluxes of potential vorticity, which requir e themselves the presence of large scale potential 
vorticity gradients, see e.g. IVallisI (|200<Th . 



Another source of fine-grained potential vorticity would be provided by the addition 
of bottom topography. This should in this case play against barotropization, since the 
topography induces potential fine-grained enstrophy in the lower layer only. In fact, an 
initially surface intensified flow may evolve towards a bottom trapped current above the 
topographic anomaly, whic h can be expla ined by the statistical mechanics arguments 



presented in this paper, see IVenailld (|2012l ). 



We now discuss how these results may apply to the mesoscale ocean, and in what 
respect they may provide interesting eddy parametrization. First, we have considered 
the relaxation of an initial condition by an inviscid flow. Real oceanic flows are forced 
dissipated. Statistical mechanics predictions should apply if a typical time scale for self- 
organization is smaller than typical time scale of forcing and dissipation. Issues of forcing 
and dissipation lead to another difficulty: what is the domain in which the flow self- 
organizes ? Strictly speaking, it should be the oceanic basin. But forcing and dissipation 
become dominant at basin scale (with, for instance, the Sverdrup balance setting the gyre 
structures). In practice, on a phenomenological basis, one could still apply the statistical 
mechanics predictions by considering an artificial domain scale with prescribed scale of 
arrest for the inverse cascade, governed say by dissipation processes, but which could not 
be predicted by statistical mechanics itself. 

We have also neglected bottom friction, because equilibrium statistical mechanics apply 
for non-dissipative systems. It is gene rally believed that bot t om friction plays a key role i n 



the vertical structure of oceanic flows lArbic fc Flierll (|2004l ): lThompson fc Youn 



ey role m 
3 (|2006h . 



On a phenomenological basis, quasi-statistical equilibria could be computed in the high 
bottom friction limit by adding a constraint of zero velocity at the bottom, for instance 
by considering ■ip(-H) = in equation ()3.14j) . 

Third, equilibrium statistical mechanics predict the final state organization of the freely 
evolving inviscid dynamics, but neither predicts the route toward this state nor the corre- 
sponding typical time scales for the convergence toward this state. Various parametriza- 
tions that relax the flow toward the equilibrium state, following a path of maximum en- 
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tropy production have been proposed, see IKazantsev et~~al. ( 19981 ): Robert fc Sommeria 
(|l992h . These parametrization satisfy basic physical constraints satisfied by the dynamics 
(PV distribution conservation and energy conservation), but the actual dynamics may 
in some cases f ollow a different path than one of maximu m entropy production. 

For instance, ISmith fc Vallis! (l200lh : iFu fc Flierll <|l980l) reported the existence of two 
time scales for the flow evolution in a freely evolving quasi-geostrophic turbulent flow with 
surface intensified stratification. On a fast time scale, the dynamics lead to a inverse cas- 
cade of vertical modes toward the first baroclinic modes, which is surface intensified, and 
observed the tendency toward barotropization at much larger time scale. The existence 
of these two time scales is a dynamical effect that can not be accounted by statistical 
mechanics. 

The results presented in this paper have highlighted the important role of the fine- 
grained enstrophy profile, for this imposes strong constraints for the eddy structures. 
On a practical point of view, from the perspective of eddy parametrizations, our result 
suggest that rather than assuming a vertical structure of eddies given by SQG modes, 
or given by combination of barotropic and first baroclinic modes, one could compute the 
eddy structure with equation (|3.14|) . assuming that K is the scale of arrest, and Zq an 
enstrophy profile that could be deduced from the resolved flow. For instance, Zq could 
be set by the structure of the most unstable mode of local baroclinic instabilities. 

Observations indicate the presence of long-lived surface intensified eddie s in the ocean 
( Chelton et oZll2007t I Johnson fc McTaggartl feoiO: IChaigneau et all 120111 ). and we here 
speculate as to how these eddies may be interpreted according to a statistical mechan- 
ical theory. First, the eddies may be interpreted as local statistical equilibrium states 
of the continuously stratified dynamics associated with a surface intensified fine-grained 
enstrophy profile, much as we considered in the present paper. Second, they could be 
interpreted as statistical equilibrium states of a 1.5 layer qu asi-geostrophic mode l assum - 
ing that only the upper layer of the ocean is active, see e.g. IVenaille fc Bouchet ( 2011a ). 
Third, they could be far from equilibrium states driven by physical mechanisms prevent- 
ing convergence towards statistical equilibrium (dynamical effects due to non-uniform 
stratification, bottom friction, large scale forcing). Understanding which hypothesis is 
relevant to explain the formation of long-lived surface intensified ocean vortices will re- 
quire further research. 
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Appendix A. Minimum enstrophy states 

The aim of this appendix is to show that the class of MRS statistical associated with 
a linear q — tp relation can be computed by solving a minimum enstrophy variational 
problem satisfie d by the coarse -grained state q. This is a dire c t consequen c e of a more 
general result of lBouched (|2008h . see also lMaida fc Wangl (^OOd k lNaso et all (|2010h . Here 
bottom topography is omitted to simplify the presentation, but taking into account this 
additional parameter is straightforward. In addition, the results obtained in this appendix 
are valid whatever the domain geometry. The only complication added by the presence 
of boundaries is that the circulation at these boundaries becomes a parameter that must 
be discussed. 
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A.l. Variational problem of the MRS statistical theory 

We consider the variational problem (|3.8p given by MRS theory. To compute its critical 
points, we introduce the Lagrange multipliers /3 t ,a(a, z),C,(x,y) associated respectively 
with the energy £ (|3.7p . with the global PV distribution V (|3.4p . and with the local 
normalization constraint J\f given by (|3.ip . Critical points are solutions of 

SS - t 8£ -(^J dz J da aSV^j - (j dxdy (SAf^j = 0, 

where first variations are taken with respect to the field p. This equation and the nor- 
malization constraint lead to 

1 . f 

exp (/3tV , c — ct (cr, z)) ,withz Q (w)= / da exp {ua — a (cr, z)) . (Al) 



The coarse-grained potential vorticity field associated with this state is 

q = J dcr ap = f a (J} t il>) , with f a (u) = — In z a ; (A 2) 

and the local fluctuations of potential vorticity are 



q>-q 2 = I daa 2 p = f a (M), (A3) 
where f' a is the derivative of f a . 



A. 2. Case of a linear q — ip relation 
Here we focus on a particular class of solutions, assuming 

f a (u) = Q(z)u, 

where the coefficient Q(z) will be expressed in term of the constraints of the problem. 
Equations (|A 2[) and (IA 3p give the q — ip relation and the fluctuations 



q = /3 t n(z)iP, n = q 2 -q 2 . 
Integrating the fluctuations on the horizontal leads to 

il= { dxdy (q*-q 2 ) =2(Z-Z cs ),with (A 4) 



v 



Z ^ = \\ T) dxdy (X dfJ a2p ) ' Zcg ^ = \ J v dxdy (X d(J ap ) ' 
Using the relations (|A 111 A 2ft . we find that 

z Q H = (2^) 1 / 2 ex P Q^ 2 ), p=-^ exp (- { ^f-Y a(a,z) = ^. 



(A 5) 



A. 3. Link with a minimum enstrophy problem 



Solution of dual, less constrained variational problem i s always a soluti on of the more 
constrained problem for some values of the constraints ( Ellis et al. 2000h . We relax the 



constraints on the PV distribution by introducing the Legendre Transform Q a of the 
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entropy associated with this constraint, taking into account the conservation of the total 
energy Eq and of the fine-grained enstrophy Zq{z) in each layer: 

G(E ,Z (z),a) = min {Q a [p] \ £ [p] = E , Z[p] = Z (z)} , (A6) 
p,N{p]=l 

Q a = I dxdy I dz da phi p + dz da a dxdy p, (A 7) 

Jv J-h J J-h is Jv 

where £ is given by (|3.7p , where the pdf field is locally normalized using (13. ip . Equations 

(|A~5|) and (TA~4|) give 

/ x ^ , CX P (- 4(Z -C) ) 
a \ a i z ) = 7777 7^7" and P = 777 

4(Z a -Z cg ) (4n(Zo - Z cg )) 1/2 

Injecting these expressions into (|A 71) yields 



Get [p] = - / dz ic [ ) + est, with = — ln(l 

2 J-h V / 



x) + 



In the limit x <C 1, we have w a;, so finding the minimizer p m of C/ Q with 

the (global) constraint on the energy and the enstrophy is equivalent of finding the 
minimizer of the functional /_„ dz Z cg [q]/Zo with the same constraints. The hypothesis 
that Z cg <C Zq is consistent with the "strong mixing limit" ((3 t ipa 1) or the low energy 
limit for which MRS equilibria are characterized by a linear relation between potential 
vorticity and streamfunction. We conclude that MRS equilibrium states characterized by 
a linear q — tp relation are solutions of the variational problem 
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